Tracing the thermal mechanism in femtosecond spin dynamics 
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We compare femtosecond pump-probe experiments in Ni and micromagnetic modelling based on 
the Landau-Lifshitz-Bloch equation coupled to a two-temperature model, revealing a predominant 
thermal ultrafast demagnetization mechanism. We show that both spin (femtosecond demagnetiza- 
tion) and electron-phonon (magnetization recovery) rates in Ni increase as a function of the laser 
pump fluence. The slowing down for high fluences arises from the increased longitudinal relaxation 
time. 
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The implementation of novel magnetic recording and 
spintronic devices requires well-funded knowledge of the 
limits of spin manipulation. Pump-probe experiments 
with powerful femtosecond lasers |l|, 0, H, U 0, 0, 0, 
Q have pushed these limits down to the femtosecond 
timescale in the past decade. These experiments have 
attracted many researchers with the aim to understand 
both fundamental mechanisms of the magnetization dy- 
namics in a strongly out-of-equilibrium regime and to 
control the magnetic properties of materials on the fem- 
tosecond timescale. Very recently the involvement of the 
spin in thermoelectric processes, spin Peltier or Seebeck 
effects, has become of strong interest Q. Therefore a vi- 
tal understanding of energy transport between the spin, 
electron and phonon system in ferromagncts is needed. 
Even for the most simple itinerant ferromagncts, such as 
Ni, the processes connecting the elementary spin scatter- 
ing process and the terahertz (THz) spin- wave generation 
have not been identified yet. However, they are the key 
to understand the macroscopic demagnetization on the 
femtosecond time scale. 

Different non-thermal mechanisms of how light could 
couple to the spin system have been put forward: the ex- 
citation of non-magnetic states mediated by the enhanced 
spin-orbit interaction (SOC) [lot . the inverse Faraday Q 
or the Barnet effect At the moment, no clear proof 
of these effects has been presented for the ultrafast mag- 
netization dynamics in Ni: (i) The change of the light 
polarization has no influence on the femtosecond demag- 
netization and estimations suggest that the amount of di- 
rect angular momentum transfer from photons to spins is 
negligible (ii) the time-dependent density functional 
theory based on the SOC mechanism at its present state 
of the art srongly underestimates the experimentally ob- 
served timescales in itinerant ferromagnets [l^l ■ 

Differenty from the latter, it has been shown that more 
applicable is a "thermal" ansatz for the description of 
the femtosecond magnetization dynamics in the itiner- 
ant ferromagnets [11,11, H, 13,[l4|. Within this descrip- 
tion, it is assumed that the excited state is a statisti- 
cal ensemble of many electronic excitations, based on 
the undisputed fact that the photons of a femtosecond 



laser, focused on a metal, pass the energy to the sub- 
systems of electrons, phonons and spins. Photons are 
absorbed by electrons close to the Fermi level leading to 
a non-equilibrium distribution that thermalizes within 
several femtoseconds. In the so-called two-temperature 
(2r) model [il,|il,[l3l, 

energy dissipation in the system 
is introduced by including rate equations, artificially sep- 
arating electrons and lattice by attributing temperatures 
to both. The stochastic spin-flip processes are driven by 
the increased electron temperature. The most probable 
candidate for a microscopic, femtosecond spin-flip mecha- 
nism is the phonon- and impurity-mediated Elliott- Yafet 
scattering The corresponding spin-flip dy- 

namics model should be able to directly relate the elec- 
tronic excitations and the magnetic precessional motion 
arising from a Landau-Lifshitz type damped precession. 
In this line a thermal formalism based on the stochas- 
tic Landau-Lifshitz-Gilbert (LLG) equation was used to 
model the sub-picosecond spin dynamics within an atom- 
istic Heisenberg model by Kazantseva et al. [l^. The 
spin relaxation rates calculated from these models have 
been found to depend on the stren gth of the coupling to 
the electron temperature bath H, 23| 
presence of THz high energy spin waves 0, . 

Different to the discrete Heisenberg model, we use 
in the foUwing a micromagnetic (macrospin) Landau- 
Lifshitz-Bloch (LLB) equation where the magnetic fluc- 
tuations are described by using the temperature and 
time-dependent parameters. We compare this thermal 
macrospin model with the measured ultrafast magneti- 
zation and reflectivity dynamics in Ni for a set of differ- 
ent excitation powers. The excellent agreement between 
experiment and modelling proves that the ultrafast de- 
magnetization in Ni has a purely thermal origin. 

Our experiment was performed on a 15 nm thick ferro- 
magnetic Ni film deposited on a Si(lOO) substrate us- 
ing electron beam evaporation in ultra high vacuum 
(5 X 10^^°mbar). A femtosecond pump-probe experi- 
ment was used to determine the reflectivity R(t) and the 
Kerr rotation Ok(t) as a function of the probe pulse delay 
T [3, 01 immediately after sample preparation. The excit- 
ing pump pulse fluence was varied from 10 to 50mJ/cm^ 



and solidifled the 
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per pulse (80 fS; A = 800 nm) calculated from the area 
determined from the half width of the Gaussian pump 
beam intensity profile (60 /xm). A static field was ap- 
plied in-plane to saturate the sample. The Kerr rota- 
tion Ok(t) is defined in the following as the asymmetric 
part of the signal, 9k~{t) = ^{9k{t, M) - 0k{t, -M)) 
changing with the field direction and mirroring predom- 
inantly the magnetization change. The symmetric part 
OK.+ ir) = ^{0k{t,M) + Oxir, -M)) mirrors the reflec- 
tivity change R{t). However the reflectivity change R{t) 
was determined separately by a differential measurement 
method of the direct sample reflectivity for p-polarized 
probe pulses. 
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FIG. 1: (a) Response after femtosecond laser excitation re- 
vealing the dynamics of the heated electrons, coherent stress 
waves and incoherent lattice excitations at various pump flu- 
ences F contributing to the transient reflectivity for a 15 nm 
Ni thin film, (b) The results of the integration of the 2T 
model for two pump fiuence values. 

The data of the delayed reflectivity R{t) for various 
pump flucnccs arc presented in Fig. [IJa) and show os- 
cillations with a main period of about 7 ps. It originates 
from the coherent excitation of a Shockwave (sound ve- 
locity vjssi = 5400 m/s). In case of a film thickness in the 



range of tNi 



opt, 800 nm 



a standing shock wave forms 



with the wave vector k = it /tm [HI. For a threshold fiu- 
ence above 30 mJ/cm^ the strong coherent stress wave ex- 
citations appear to have an incoherent counterpart (lat- 
tice heating and thermal expansion) owing an effect on 
the reflectivity. 

Neglecting the contribution of the incoherent stress 
wave excitations, the reflectivity change in the first ap- 
proximation is defined by the electron temperature Te 
[Iq . The electron temperature Te is coupled to the lat- 
tice temperature T/ within the 2T model |l5| in the form 
of two differential equations: 



= -G,i{T,,-Ti)+P{t) 

at 

at Tth 



Here Ce and Ci are the specific heats of the electrons and 
the lattice, Gei is an electron-phonon coupling constant 
which gives the rate of the energy exchange between the 
electrons and the lattice [l3| and Tth is the heat diffu- 
sion time. The source term P{t) is a function which de- 
scribes the laser power density absorbed in the material. 
We will assume a Gaussian laser pulse of 80 fs duration. 
The electron specific heat should strongly depend on Te 
due to variation of the density of states of Ni around 
the Fermi level with temperature. Here we use a sim- 
plification for this dependence and assume Ce = Te^e, 
7e = 3 X 10'^ Jm^'^K^^. Within this approximation the 
analysis of the data in Fig. [IJa), using the method of 
Ref.Ill, gives Gei « 10 x lO^Wm-^K"! at different 
levels of laser excitation, consistent with similar values 
reported in literature for Ni [2^ . 
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FIG. 2: Parameters defining the slowing down of the magne- 
tization rates, (a) The ratio between the maximum electron 
and lattice temperatures 1^^°-^ as obtained fitting the 

2T model to the refiectivity data in Fig. [ija) as a function 
of pump fiuence. (b) The longitudinal relaxation time for 
Ni as a function of temperature evaluated from the MFA. 
(c) Schematics of the thermal spin excitations followed by 
the femtosecond laser pulse (left) and the thermal macrospin 
model (right). 

The lattice specific heat at low pump fluences is taken 
as Ci = C(300 K) - Ce(300 K) = 3.1 x 10^ Jm-^K"!, 
where C(300 K) = 4 x 10^ Jm~^K"^ is the experimen- 
tal total specific heat of Ni The fitting of the 2T 
model to the reflectivity data from Fig. [TJ^a) shows that 
the energy, put into the lattice, increases as a function 
of pump fluency, showing a decrease of the Tmax/^maa: 
ratio (see Fig. Wp)- This decrease may contain different 
contributions such as (i) nonlinear temperature depen- 
dence of the specific heat Ce, (ii) additional energy dis- 
sipation mechanisms, e.g. excitation of elastic waves for 
T > 6Debye (475 K in Ni) (iii) dependence of the refiec- 



(1) 



tivity data on T/ at high pump fluencies [2l|. The 2T 
model was adjusted to capture the behavior observed by 
decreasing the lattice specific heat from the value above 
to Ci = 1.0 X 10^ Jm~^K~^ at high pump fiuence. Fi- 
nally, the heat diffusion time of Tth ~ 50 ps was obtained 



3 



(a) 



1.0 



I 



E 0.9 



< 



< 0.7 



F=10mJ/cm 




(b) 



300' 
■ 200' 
= 100' 



■ experiment 




^ LLB model 




. ♦ 8 1 













Ax [ps] 



20 40 60 

Fluence [mJ/cm ] 



FIG. 3: (a) Absolute values of the measured dynamics of the 
Kerr angle rotation (solid lines) and the modeled magnetiza- 
tion dynamics (symbols) for various laser pump fluences F. 
(Right) The magnetization dynamics rates as a function of 
pump fluence: (b) demagnetization time tm and (c) mag- 
netization recovery time te as function of the pump fluence 
are plotted. Experimental (squares) and modeling (circles) 
values are given for a direct comparison. Lines are a guide 
to the eyes. Note that te differs from the electron-phonon 
equilibration Tep by its magnetic contribution. 



by fitting the reflectivity data in Fig. [l]for all pump flu- 
ences. Figure presents the result of the 2T model 
integration (Eqs. ([1])). Note that for the highest pump 
fluence the electron temperature exceeds the Curie tem- 
perature Tc^Ni while the lattice temperature stays well 
below that value. 

The magnetization dynamics model is based on the 
LLB micromagnetic equation [2^ . This macroscopic 
equation has been derived for thermodynamically aver- 
aged spin polarization m by D. Garanin within mean 
field approximation (MFA) from the classical Fokker- 
Planck equation for atomistic spins interacting with a 
heat bath and from the corresponding density-matrix in 
the quantum case (25j . The macrospin LLB equation 
has been shown to be a valid micromagnetic equation 
at all temperatures, even above Tc [1^. The advantage 
of the present model is the use of only one thermody- 
namically consistent macrospin equation, as compared 
to large-scale calculations based on the atomistic Heisen- 
berg model We write the LLB equation as follows: 



rh = 7[m x Hoff] + — ^^-(m • (Hcff))m 



7a_L 



m X [m X (Hoff)]], 



(2) 



where 7 is the gyromagnetic ratio, a|| and a±_ are dimen- 
sionless longitudinal and transverse damping parameters 
given by a\i = 2XT/3Tc, a_L A [1 - T/STc] for T < Tc 
and a± = a\\ for T > Tc- Here A is the parameter de- 
scribing the coupling of the spins to the heat bath, which 
is assumed to be the electron system for the present case. 



as suggested in Rcf. The effective field HoS is given 

by 



H + H, 



1 

2X1 



XII V 



1 



3Tc 



5(T-Tc) 



T<Tc 
m T>Tr 



(3) 

In our model the temperature-dependent equilibrium 
magnetization 7Tie and longitudinal susceptibility X|| were 
determined in MFA [gj], adjusting the Curie temperature 
value to the experimental one Tc = 631 K. In a more rig- 
orous multiscale description [13] , they could be evaluated 
directly from the spin Hamiltonian, parameterized from 
the electronic structure calculations, or considered as an 
experimental input. We use the following parameters for 
Ni: The saturation magnetization AI^ = 500emu/cm"^ 
(at T = OK), the applied field H = IbOOe^e, the 
damping parameter a± = 0.04 at T = 300 K [7| (from 
which the intrinsic coupling parameter A = 0.045 is esti- 
mated). The easy plane anisotropy, defined by magneto- 
static interactions for the thin film geometry, is given by 
Ha = "i^T^Ms {m^ex + myBy). The magneto-crystalline 
anisotropy is neglected. 

One of the immediate consequences of the thermal dy- 
namics based on the LLB equation is the presence of 
the longitudinal relaxation time, i.e. the rate at which 
the spin system disorders at a given temperature. The 
longitudinal relaxation time mirrors the appearance of 
short wavelength spin exciations at THz frequency in 
the thermal microscopic model. In the absence of the 
anisotropy the longitudinal relaxation time behaves as 
T||(T) = x\\{T,'H)/-fai\iT) where Xlll^'^Hlis the longi- 
tudinal susceptibility at applied field H [26|. The de- 
pendence of the longitudinal relaxation time on the tem- 
perature for Ni is presented in Fig. [2jb). It strongly 
increases at Tc^Ni, due to a faster divergence of the lon- 
gitudinal susceptibility as compared to the variation of 
the longitudinal damping; a general critical phenomenon 
known as the critical slowing down 27| . It is in the order 
of 0.1 ps for low temperatures and increases by more than 
one order of magnitude towards Tc- 

Therefrom one concludes that during the laser-induced 
demagnetization, the electron temperature is changed 
on the timescale faster than the longitudinal relaxation 
time. Consequently, the spin system cannot follow the 
electron temperature. The magnetic response is delayed 
by the slowing down of the longitudinal relaxation being 
inversely proportional to the longitudinal damping con- 
stant a||. Furthermore, a||(T) is related to the Gilbert 
damping a±{T) via the coupling parameter A. There- 
from one derives a temperature-dependent correction to 
the model of B. Koopmans et at relating the fem- 
tosecond demagnetization time to the Gilbert damping 
parameter by tm ~ l/ci_L- 

For a quantitative comparison the results of the in- 
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tegration of the system of Eqs. ([T|) and ([2|) are pre- 
sented in Fig. [3] (a) in comparison to the experimen- 
tal time-resolved Kerr angle rotation data. Both show 
the slowing down of the time scales as the pump flu- 
ence is increased from lOmJ/cm^ to 50mJ/cm^. Fig. 
[3I^a) shows a very good agreement between theory and 
experiment, giving the absolute value of the demagneti- 
zation as 3% for the lowest pump fluence to around 40% 
for the highest. To accurately determine the relaxation 
times, the slopes are analyzed using the analytic solu- 
tion for the three-temperature-modcl (including the spin 
temperature) given by Dalla Longa et al. 
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Here *G{t) represents the convolution product with the 
Gaussian laser pulse profile, Q{t) is the step function, Aj 
are the fitting constants and tm and te define the scale of 
the magnetization loss and magnetization recovery time, 
respectively. The results are presented in Fig. and 
c) . We observe a slowing of both magnetization rates tm 
and Te as a function of laser pump fluence, following the 
slowing down of the longitudinal relaxation time. As a 
result, the electron-phonon equilibration time extracted 
from the reflectivity (rep) and from the magnetization 
recovery time (te) have different magnitudes. This ob- 
servation was also reported for Fe 0]. However, for the 
slowing down of the magnetization recovery, the longitu- 
dinal relaxation slowing down is not sufficient itself. In 
Ref. [isj the experimentally observed slowing down of te 
in CoPt was attributed to the loss of the magnetization 
correlations for high pump fluences. In our experiment 
the combined measurements of static and dynamic Kerr 
effect have allowed us to determine that the demagneti- 
zation for the highest pump fluency was around 40%, not 
sufficient to produce the magnetization correlation loss. 
Our detailed analysis of the reflectivity data reveals that 
at high pump fluences, the relative energy loss is faster in 
the electron system than in lattice, mirrored by the ratio 
between maximum electron and lattice temperatures pre- 
sented in Fig. IHa): The cooling dynamics of the lattice 
and electrons, together with the longitudinal relaxation, 
leads to the observed slowing of the magnetization recov- 
ery. To summarize the result of the direct comparison, 
the thermal macrospin model, based on the concept of a 
longitudinal relaxation time in combination with the 2T 
model, fully describes the absolute variation and charac- 
teristic time scales of the femtosecond demagnetization 
experiment in Ni. 

In conclusion, the agreement between the thermal 
macrospin model and the experimental data reveals 
the dominant character of the thermal demagnetization 



mechanism in Ni mediated by hot electrons. This thermal 
magnetization dynamics model has been used to quan- 
titatively describe the absolute values and characteris- 
tic timescales determined in the ultrafast magnetization 
and reflectivity data with increasing the excitation flu- 
ence. As a result the theoretical model allows us to iden- 
tify the longitudinal relaxation as a limiting factor for 
the slowing down of the demagnetization rates, a general 
phenomenon which should be present in all ferromagnetic 
materials. Similarly, other 3d-metals such as Fe, where 
the high temperature dynamics of the itinerant electrons 
is suitable for the heat bath description, can be described. 
However, at the same time the longitudinal relaxation 
rate is also dependent on the coupling strength between 
electron and spin. The latter means that the magnitude 
of the demagnetization rate strongly depends on the mi- 
croscopic electron, lattice and spin coupling mechanisms. 
This suggests that femtosecond pump-probe experiments 
are a valuable tool to characterize spin based caloric ef- 
fects of a particular material. 
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